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PREFACE 
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A(t)  *  age  at  time  t  along  the  ray 

A'(xj_)  *  area  distribution  of  slender  body 

C  *  pressure  coefficient 

F  *  F-function  for  aircraft  signatures 

FA  3  area  component  of  F-function 

Fq  =  lift  component  of  F-function 

F  ,  ,  F  „  =  nondimensional ized  area  and  lift  contributions 

A1  B 1 

Ff  =  F-function  conversion  factor 

F ^  =  input  F-function 

j  =  Jacobian,  used  to  define  raytube  area 

Ks  =  aircraft  shape  factor 
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Lp  =  a  length  used  to  nondimensionalize  the  F-function 

Lce  =  C -weigh ted  sound  exposure  level  (CSEW,  dB) 

Lpk  =  maximum  flat-weighted  sound  pressure  level,  dB 

M  3  aircraft  Mach  number 

0  *  aircraft  heading 

p  =  atmospheric  pressure 

PA  =  atmospheric  pressure  at  the  aircraft 

p(<j>,i|,)  =  pressure  far  field  signature 

R  =  relative  radius  of  curvature 

S  distance  from  the  caustic  along  the  ray 

»  area  distribution  of  equivalent  body  of  resolution 
V  =  aircraft  speed 

VA  =  average  velocity  vector 

W  =  aircraft  weight 

X*  =  distance  to  focal  zone  boundary,  tangential  to  the  focus 

y  =  distance  normal  to  the  focus 

Y*  =*  distance  to  focal  zone  boundary,  normal  to  the  focus 

c  *  local  speed  of  sound  at  a  gi,7en  altitude 

cA  =  speed  of  sound  at  aircraft 

1  =  aircraft  characteristic  length,  typically  fuselage 

tough  (Section  4) 
nL  =  lift  load  factor 


Subscr  ipt 
A  =  aircraft 

o  -  initial  value  at  time  of  emission  of  a  ray  from  aircraft 
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ny > 
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x^A) 

z(tA) 


4u4  V 


<  R 
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=  unit  normal  components  in  x,  y,  and  z  direetionb 
*  vector  components  of  normals  to  phase  surfaces  of  wave 
numbers  (relative  to  airborne  reference  frame) 

=  atmospheric  pressure  at  aircraft  altitude  (Section  4) 

=  .5PAy  My2  =  dynamic  pressure 
=  distance  from  aircraft  flight  path 
=  aircraft  speed  at  time  t 
=>  time  along  ray 

=  wind  components  at  a  given  altitude 
=  coordinate  system  near  aircraft  aligned  with  ray 
=  ground  coordinate  system  (east,  north,  and  height 
above  sea  level,  respectively) 
tA)  / 

=  aircraft  location  at  time  tA 
=  (Y  +1  )/2 

=  local  perturbation  velocity  potential  function 


=  ratio  of  specific  heats  for  air 
=  airspeed  component  in  horizontal  plane 
=  peak  overpressure 
*  aircraft  airspeed 

=  perturbation  pressure,  incremental  pressure  to  the 
sonic  boom 

=  incremental  change  in  ray  location  =  800  ft 
=  velocity  difference  of  wind  components  at  aircraft 
altitude 

3  caustic  curvature  vector 
=  ray  curveture  vector 

=  relative  curvature  vector  between  ray  and  caustic 


surface 

*  Mach  cone  angle,  sin~l(  ) 
=  position  behind  the  nose 


Subscript 
A  =  aircraft 

o  =  initial  value  at  time  of  emission  of  a  ray  from  aircraft 
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p  *  atmospheric  density 

$  =  azimuth  angle  of  ray  from  vertical  plane  or  ray  bank 

angle 

<p  inc  =  incremental  ray  bank  angle  per  unit  distance 

f  =  phase  (identified  with  position  ) 

u  =  frequency  (scaled  by  A/C  length)  equal  to  the  airspeed 

in  the  airborne  reference  frame 


Subscript 
A  *  aircraft 

o  =  initial  value  at  time  of  emission  of  a  ray  from  aircraft 
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1.0  INTRODUCTION 


The  BOOMAP2  and  MOAOPS  computer  programs  are  utilized  to 
analyze  noise  from  supersonic  aircraft  operations  in  Military 
Operating  Areas  (MOA's).  The  two  programs  are  designed  to  extract 
and  analyze  information  from  the  Tactical  Air  Crew  Combat  Training 
System/Air  Combat  Maneuvering  Intrumentation  (TACTS/ACMI)  manu¬ 
factured  by  the  Cubic  Corporation.  The  TACTS/ACMI  system  digi¬ 
tizes  various  positional  and  performance  parameters  of  the  air¬ 
craft  in  a  Military  Operating  Area  at  frequent  intervals  for  later 
replay  in  graphic  or  tabular  form  during  air  crew  briefings. 

The  MOAOPS  program  (Ref.  1)  extracts  information  from  a 
TACTS/ACMI  mission  standard  data  tape  and  compiles  a  computer 
library  of  information  concerning  the  supersonic  operations.  The 
B00MAP2  program  utilizes  the  library  produced  by  the  MOAOPS 
program.  The  program  calculates  various  statistics  on  the  super¬ 
sonic  operations.  It  also  calculates  expected  sonic  boom  levels 
on  the  ground  based  on  the  extracted  information.  Both  programs 
are  written  in  FORTRAN  77  and  operate  in  batch  mode  on  Control 
Data  Corporation  (CDC)  CYBER  170  Series  machines. 


The  BOOMAP2  program  is  capable  of  predicting  the  noJse  levels 
or  overpressures  on  the  ground  resulting  from  either  carpet  or 
focus  sonic  booms  resulting  from  air  combat  maneuvering  training 

C  "1  J  A,  _  rT»L  •  L  ^  W  U  r*  Art  O  1  Q  f  C  Qr  ! 

i.  .L  iQllUd*  Ail©  uuuyuu  w  -  • 


(a)  Various  statistical  summaries. 

(b)  Flight  track  information. 
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(c)  A  computer  library  of  predicted  overpressures  on  the  ground 
for  each  flight  analyzed. 

(d)  "Scratch  pad"  plots  showing  maximum  overpressures  when 
focused  sonic  booms  occur. 

(e)  Calculated  noise  levels  at  a  grid  of  ground  positions 
{  100  x  100  matrix ? . 

The  output  of  the  BOOMMAP2  program  is  designed  to  be 
compatible  with  GPCP  (General  Purpose  Contouring  Program)  (Ref. 

2).  Through  the  use  of  the  GPCP  program,  the  BOOMAP2  output  can 
be  displayed  in  terms  of:  (a)  a  map  showing  the  flight  track 
segments  where  the  aircraft  was  supersonic?  (b)  displays  of  the 
calculated  sonic  boom  "noise"  in  terms  of  several  metrics;  and  (c) 
a  map  showing  the  location  of  focused  sonic  boom3. 

The  BOOMAP2  computer  model  utilizes  a  sophisticated  acoustic 
ray  theory  model  for  predicting  the  overpressures  and  noise  levels 
on  the  ground,  which  is  based  upon  the  TRAPS  program  developed  by 
Dr.  Albion  Taylor  (Ref.  3).  The  version  of  tho  TRAPS  program  in 
BOOMAP2  incorporates  several  corrections  and  changes  over  the 
original  program.  The  BOOMAP2  program  also  incorporates  several 
additions  to  the  original  TRAPS  program  to  permit  estimation  of 
overpressures  at  focus  locations. 

The  BOOMAP2  program  replaces  the  BOOMMAP  program  developed 
earlier  (Ref.  1).  The  major  difference  betwee  programs  is  that 
the  original  BOOMMAP  program  utilized  calculation  procedures  for 
estimating  booms  based  on  the  simplified  sonic  boom  prediction 
model  methods  developed  by  Henry  Carlson  (Ref.  4).  These  calcula¬ 
tions  assumed  steady  state  flight  and  were  not  intended  to  provide 
accurate  predictions  for  focus  booms.  The  current  program  pro¬ 
vides  much  more  accurate  calculation  of  overpressures  at  and  near 
focus . 
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This  report  is  one  of  three  documents  describing  the  B00MAP2 
program.  The  other  documents  consist  of  a  user's  guide  (Ref.  5) 
and  a  maintenance  manual  (Ref.  6). 

Section  2  of  this  report  presents  an  overview  of  both  BOOMAP2 
and  the  MOAOPS  programs.  Section  3  provides  a  more  complete  tech¬ 
nical  description  of  the  sonic  boom  propagation  code  incorporated 
in  B00MAP2.  Section  4  describes  the  aircraft  F-functions  incor¬ 
porated  in  the  program.  Section  5  summarizes  the  rationale  for 
the  selection  of  the  TRAPS  program  as  a  basis  for  model 
calculations.  Results  from  this  program  are  compared  to  those 
from  other  ray  tracing  programs  and  field  measurements  in  Section 
6.  Section  7  provides  recommendations  for  future  study.  Appendix 
A  provides  some  examples  of  program  output. 
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2.0  PROGRAM  OVERVIEW 


The  major  purpose  of  the  BOOMAP2  and  the  accompanying  MOAOPS 
program  is  to  extract  and  analyze  information  from  the  Tactical 
Air  Crew  Combat  Training  System/Air  Combat  Maneuvering  Instrumen¬ 
tation  (TACTS/ACMI )  system  installed  at  various  combat  training 
military  operating  areas  in  order  to  predict  the  location  and  mag¬ 
nitude  of  sonic  boom  overpressures  on  the  ground  in  the  vicinity 
of  supersonic  flights.* 


Real  time  flight  information  is  transmitted  to  the  TACTS/ACMI 
systems  on  ground.  Among  the  data  is  real  time  information  on 
aircraft  position,  velocity  and  acceleration,  updated  at  intervals 
of  100  to  200  milliseconds.  The  MOAOPS  program  extracts  this  data 
for  the  sonic  boom  analysis  from  the  tapes  at  approximately  1.5 
second  intervals  in  order  to  minimize  both  the  time  taken  to  read 
the  tapes  and  the  quantity  of  information  to  be  stored. 

The  MOAOPS  program  is  in  two  parts:  a  data  extraction  program 
EXTRCT ,  and  an  index  deletion  and  modification  program  DELETE. 

The  data  extraction  program  reads  the  ACMI  tapes,  extracting  rele¬ 
vant  information  and  appending  this  information  to  either  a  new  or 
existing  data  base  (library).  The  library  file  accumulates  the 
information  from  all  the  mission  tapes  analyzed.  This  library 
file  is  indexed  so  that  a  particular  mission,  aircraft  type,  etc. 
can  be  accessed  by  the  sonic  boom  analysis  programs. 


*  In  this  report,  overpressure  will  typically  mean  the  "magnitude" 
of  the  sonic  boom  at  a  given  point  expressed  in  terms  of  the  maxi¬ 
mum  overpressure  in  pounds  per  square  foot  (psf)  or  in  terms  of 
the  overall  sound  pressure  level  (OASPL)  in  dB,  or  in  terms  of  the 
C-weighted  sound  exposure  level  (CSEL)  in  dB.  Program  options 
allow  a  choice  of  either  of  these  three  metrics  for  the  contour 
presentations . 


Functional  Relationship  Between  Elements  of  BOOMAP2  Computer  Program 


The  B00MAP2  data  analysis  program  accesses  the  MOAOPS  library 
tapes  as  selected  by  the  user.  The  data  analysis  program  produces 
statistical  and  graphical  output  describing  the  aircraft  positions 
parameters  as  various  measures  of  predicted  boom  strength.  The 
BOOMAF2  program  produces  tabular  output  of  various  statistics  that 
is  sent  directly  to  a  line  printer.  The  overpressures  predicted 
on  the  ground  by  the  acoustic  ray  theory  model  are  output  to  the 
printer  and  also  stored  as  a  computer  library  for  future  access. 

In  addition,  for  those  situations  where  focused  sonic  booms  are 
produced,  individual  plots  of  the  maximum  overpressures  together 
with  other  technical  information  are  produced  in  form  of  a 
"scratch  pad".  These  "scratch  pads"  can  be  plotted  for  each 
situation  in  which  focused  booms  occur. 

When  a  mission  is  selected  from  the  HOAOPS  library  and  used  as 
a  B00MAP2  computer  program,  the  rays  traced  by  B00MAP2  are  saved 
in  a  RAYS  library.  If  that  same  mission  is  selected  at  a  future 
time,  the  necessary  ray  information  is  recalled  from  the  library, 
thereby  saving  substantial  computer  time. 

To  produce  graphic  output,  B00MAP2  creates  a  file  which  is 
compatible  with  California  Computer  Products  (CALCOMP)  General 
Purpose  Contouring  Program  (GPCP-II).  GPCP-II  reads  this  file  and 
generates  the  necessary  plotter  directors  to  produce  hard  copy 
graphic  output. 

The  user  controls  the  data  base  subset  to  be  extracted  from 
the  MOAOPS  library  through  the  U3e  of  an  input  data  file.  Through 
this  file,  the  user  specifies:  a)  the  name(s)  of  the  MOA  ranges 
to  be  considered;  b)  mission  names  or  dates?  c)  bounding  times  of 
day;  and  d)  aircraft  types  (specific  tail  numbers  optional). 

Users  also  specify  the  desired  output  products.  These 
include : 
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1.  A  statistical  summary  of  position,  speed,  and  boom  strength 
variables.  This  summary  includes  distribution  functions  of 
range  x-coordinates  and  y-coordinates ,  and  the  aircraft 
z-coordinate  (height  above  the  range),  all  in  feet.  It  also 
includes  a  distribution  function  of  effective  height  (he). 
Distribution  functions  of  Mach  number,  cutoff  Mach  number,  and 
effective  Mach  number  are  also  presented.  Estimated  boom 
strength  distribution  functions  include  peak  overpressure  (in 
pounds  per  square  foot),  the  peak  overpressure  (in  dB,  re:  20 
microPascals ) ,  the  C-weighted  sound  exposure  level  (in  dB), 
and  the  A-weighted  sound  exposure  level  (in  dB ) .  Also 
included  are  root  mean  square  values  for  effective  height, 

Mach  number,  effective  Mach  number,  and  cutoff  Mach  number. 

2.  A  flight  track  map  depicting  ground  projections  of  flight 
paths  during  supersonic  activity. 

3.  A  flight  track  map  depicting  ground  projections  of  flight 
paths  during  sonic  boom  producing  activity. 

4.  A  noise  contour  map  of  average  C-weighted  sound  exposure 
levels  ( CSEL ) . 

5.  A  noise  contour  map  of  C-weighted  day-night  average  levels 
(CLDN ) .  This  requires  input  of  the  reference  number  of 
daytime  operations  which  is  used  to  convert  CSEL  to  CLDN. 

6.  A  noise  contour  map  of  average  peak  overpressures  in  pounds 
per  square  foot,  or  OASPL, 

7.  A  map  showing  the  geographic  location  of  maximum  overpressures 
due  to  focused  sonic  booms. 

Examples  of  the  BOOMAP2  program  output  are  shown  in  Appendix  A. 
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3.0  TECHNICAL  OVERVIEW  OF  THE  SONIC  BOOM  PROPAGATION  CODE 


This  report  provides  a  technical  overview  of  the  sonic  boom 
propagation  code  incorporated  in  BOOMAP2.  The  code  is  based  upon 
the  TRAPS  program  developed  by  Dr.  Albion  Taylor  (Ref.  3).* 
Modifications  to  the  code  consisted  of 

a)  eliminating  some  programming  errors  that  were  discovered  in 
the  course  of  the  project; 

b)  augmenting  the  program  with  a  technique  for  estimating  the 
sonic  boom  signatures  at  a  simple  focus  using  a  similitude 
developed  by  Gill  and  Seebass  (Ref.  7)  and  originally  implemen 
ted  by  Plotkin  (Ref.  8).  Portions  of  the  code  from  the  FOBOOM 
program  developed  by  Dr.  Kenneth  P.lotkin  have  been  used; 

c)  constraining  the  propagation  to  consider  only  the  portion  of 
the  sonic  boom  wave  front  originating  beneath  an  aircraft  and 
the  propagation  through  a  standard  still  (no  wind)  atmosphere; 

d)  developing  a  driver  for  selection  of  the  portion  of  the  sonic 
boom  footprint  to  be  traced  and  ground  signatures  to  be  saved. 

A  body  moving  through  the  atmosphere  at  supersonic  speeds  will 
continuously  generate  a  system  of  shock  waves  in  its  wake.  Under 
appropriate  atmospheric  conditions  these  sonic  boom  shock  waves 
will  produce  a  disturbance  at  the  ground.  Computer  programs  for 
estimating  the  ground  level  sonic  boom  must  include  the  following 
elements ; 

•  A  method  for  generating  (describing)  a  trajectory  for  the 
supersonic  craft 


*  Much  of  this  section  is  adopted  directly  from  Ref.  3. 
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•  A  method  for  characterizing  the  system  of  shock  waves 
generated  about  the  craft 

•  A  description  of  those  atmospheric  parameters  that  affect 
the  propagation  of  the  sonic  boom  shock  waves 

•  A  method  for  characterizing  the  path  of  the  shock  wave 
propagation  through  the  atmosphere 

•  A  method  for  evaluating  the  effect  of  the  atmosphere  on 
the  magnitude  and  shape  of  the  sonic  boom  waves 


•  A  method  for  accounting  for  the  reflection/attenuation  of 
the  sonic  boom  at  the  ground 

•  A  technique  for  assembling  and  presenting  the  results  in 
a  meaningful  fashion 


3,1  The  Aircraft  Flight  Path 


Both  the  original  TRAPS  program  and  BOOMAP2  use  recorded 
flight  path  information  to  characterize  the  supersonic  aircraft 
trajectories.  In  TRAPS,  the  displacements  of  the  aircraft  are 
used  to  fit  a  cubic  spline  from  which  the  acceleration  vectors  are 
calculated.  These  are  smoothed  and  an  inverse  spline  is  used  to 
calculate  the  smoothed  aircraft  locations. 


B00MAP2  uses  a  standard  cubic  spline  fit  of  the  velocities  to 
calculate  the  acceleration  vectors.  The  accelerations  are 
smoothed  and  quadratic  coefficients  are  calculated  using  a 
weighted  linear  least  squares  method  which  are  then  used  to 
interpret  the  aircraft  flight  data  at  any  specific  time  between 
the  input  aircraft  track  times. 
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The  flight  path  information  serves  three  functions  in  the 

sonic  boom  propagation  analysis: 

1.  It  provides  an  initial  location  from  which  the  sonic  boom 
shock  wave  is  propagated.  At  each  time  point  along  the  flight 
path  the  sonic  boom  shock  wave  system  originates  at  the  lead¬ 
ing  edge  of  the  aircraft;  other  aircraft  features  generate  the 
detail  of  the  waveform  at  the  location  of  these  features.  The 
length  of  the  aircraft  (and  hence,  the  initial  length  of  the 
system  of  shock  waves)  is  much  shorter  than  the  distance  from 
the  aircraft  to  the  ground.  Thus,  the  path  of  the  shock  wave 
system  may  be  characterized  by  tracing  the  path  of  the  leading 
shock . 

2.  The  leading  shock  wave  (near  the  aircraft)  will  be  a  conical 

wave  (called  the  Mach  one)  with  the  axis  of  symmetry  along  the 

•> 

aircraft  velocity  vector,  VA  (Figure  2).  The  aircraft  Mach 
number,  M  (M  -  j  V^j  /cA,  cA  is  the  speed  of  sound  at  the 
aircraft),  is  related  to  the  apex  half-angle,  u  ,  by  1/M  <* 
sin  .  The  initial  direction  of  propagation  of  the  wave  is  at 
ninety  degrees  to  the  surface  of  the  Mach  cone?  the  collection 
of  the  initial  directions  of  propagation  forms  a  cone  called 
the  ray  cone.  Thus,  the  aircraft  velocity  vector  determines 
the  initial  direction  of  propagation  of  the  shock  wave. 

3.  The  ray  cones  generated  by  an  aircraft  accelerating  along  its 
velocity  vector  will  have  progressively  larger  apex  angles 
(Figure  3).  Thus,  corresponding  portions  of  the  wave  front 
generated  at  short  time  intervals  from  each  other  will  tend  to 
constructively  interfere  with  each  other  at  some  distance  from 
the  aircraft.  (The  distance  from  the  aircraft  at  which  this 
occurs  will  depend  on  the  rate  of  acceleration  and  the  manner 
in  which  the  atmosphere  modifies  the  propagation  path  from  a 
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FIGURE  3.  SONIC  BOOM  FOCUS  DUE  TO  ACCELERATION. 
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straight  line.)  Analogously,  an  aircraft  enqaged  in  a 
constant  speed  turn  will,  over  successive  time  periods,  gener¬ 
ate  ray  cones  whose  axes  are  misaligned  (Figure  4).  The  mis¬ 
alignment  of  the  ray  cone  axes  causes  the  ray  cones  to  be 
closer  to  each  other  on  one  side  of  the  cone  and  further  from 
each  other  on  the  opposite  side  of  the  cone.  This  effect  is 
maximal  in  the  plane  of  the  turn.  On  the  side  in  which  the 
ray  cones  are  tilted  toward  each  other  there  is  an 
intensification  of  the  signal,  while  on  the  opposite  side  the 
sonic  boom  levels  are  diminished.  Since  all  maneuvers  can  be 
described  as  combinations  of  these  basic  maneuvers,  the  effect 
of  a  maneuver  will  be  some  local  intensification 


(rarif ication )  of  the  sonic  boom.  The  magnitude  of  the 
aircraft  acceleration  and  jerk  will  affect  the  location  of 
this  enhanced  (subdued)  sonic  boom  and  the  degree  of 
amplification  or  diminishment. 


3 . 2  The  Near  Field  Signature 


In  order  to  extrapolate  the  sonic  boom  overpressure  signatures 
to  the  ground,  a  sonic  boom  propagation  code  must  be  provided  a 
description  of  the  disturbance  of  the  atmosphere  generated  by  the 
supersonic  aircraft.  This  disturbance  may  be  described  as  an 
overpressure  waveform  that  an  observer  near  the  aircraft  would 
measure  as  the  aircraft  flies  by,  or  equivalently,  in  terms  of  a 
theoretically  derived  F-function,  defining  the  flow  near  the 
aircraft.  Both  the  TRAPS  program  and  B00MAP2  computer  program  use 
this  latter  approach. 


This  approach  develops  the  initial  acoustic  siqnal  from  the 
aircraft  geometry  and  its  lift  distribution  which  are  used  to 
develop  a  velocity  perturbation  potential  function.  (A  velocity 
potential  is  a  function  whose  gradient  describes  the  velocity 
induced  in  the  fluid.)  After  a  number  of  simplifications,  the 
potential  function  can  be  expressed  in  terms  of  an  "equivalent 


■15- 


'>/>'  w^TSi ■fr'nv 

—www  mKHmmmmmmummiamMmMmMMmmtiam 


aaaa  jam  ^ 


Roy  Cone 
Generated 


kl 


area"  distribution  S'  (x^,  $  ) .  The  equivalent  area  distribution 
is  the  sum  of  two  terms.  The  first  is  the  derivative,  A' (xj.), 
taken  in  the  direction  of  the  aircraft  velocity  vector  of  the 
aircraft  cross-sectional  area  cut  by  a  plane  oriented  at  the  Mach 
angle,  y . 


The  second  is  given  by  — Mxj.)  where 

pV2 


8  ■  /M2-1 


p  =  atmospheric  density 


V  =  aircraft  speed 


£(X].)  »  the  rate  of  change  in  the  xi  direction  nf  the 

components  of  the  combined  effect  of  the  lift 
and  the  side  forces  in  the  negative  $  direction. 


The  resulting  asymptotic  expression  for  the  potential  function  is 

,  rx~Br 

‘{’(x-Br,  r,  $}  =  -  S'  (x^Jdx^ 

2  yt  /i  Br  -<®  , — 5 - 

H  /x-Br-x, 

where  r  is  the  distance  from  the  aircraft  flight  path. 
The  perturbation  pressure,  p,  is  given  by  -pv2  x.  On 
differentiating  the  potential  function  this  gives 
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Alternatively,  the  relationship  for  the  perturbation  pressure 
may  be  written  in  the  form 

PCR  'r 

where 

F  »  F.  F - 
i  f 

_ 

F  -  *  - 

f  m 

This  F-function  can  be  further  decomposed  into  area  and  lift 
components  FA  and  Fq  as  follows 


F  =  Ff(FA  +  Fq  cos  ) 


^  * 


Fai (Lf) 
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.5  S  (n_  W/  ) 

fb  -  fbi  — TiA 


(Lp) 


where 


nL  -  li;  t  load  factor 
W  =  aircraft  weight 
q  =  0.5Pay  m2  »  dynamic  pressure 
PA  =  atmospheric  pressure  at  the  aircraft 
Y  =  ratio  of  specific  heats  for  air 
Lp  s  a  length  used  to  nondimensionalize  the 
F-function 

=  Nondimensionalized  area  and  lift 
contributions » 


Note  that  an  unnormalized  F-function,  FA,  may  be  related  to  a 
measured  nearfield  signature  close  to  the  aircraft  by  the 
relationship 

Ap  ■  Y?aM2  Fa 
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The  TRAPS  program  input  consists  of  the  component  F-functions 
FaI  and  FbI*  In  BOOMAP2,  the  TRAPS  code  to  handle  the 
component  F-functions  has  been  retained  intact,  but  is  presently 
unused.  The  BOOMAP2  code  employs  simplified  F-functions  based 
upon  a  procedure  developed  by  Carlson  (Ref.  4).  The  development 
of  the  simplified  F-function  is  described  in  Section  4  of  this 
report. 

3.3  The  Atmospheric  Description 

While  the  propagation  of  sonic  boom  through  the  atmosphere  is 
affected  by  the  detailed  characterization  of  the  atmosphere  (tem¬ 
perature,  pressure,  winds,  and  chemical  composition),  as  a  practi¬ 
cal  matter  useful  results  may  be  obtained  by  employing  simplified 
descriptions.  A  common  approach  used  is  to  describe  these  atmos¬ 
pheric  properties  as  being  horizontally  stratified  and  temporally 
constant.  This  has  the  practical  consequences  of  making  the  prop¬ 
agation  analysis  significantly  more  tractable  and  creating  mathe¬ 
matical  models  consistent  with  the  meteorological  data  available 
to  use  with  them  while  sacrificing  only  the  ability  to  model  what 
i re  normally  rather  localized  effects  (focusing  or  attenuating 
sonic  boom  waveforms  by  small  scale  anomalies).  This  omission  is 
likely  to  be  important  only  for  those  atmospheric  variations  in 
the  immediate  vicinity  of  either  the  aircraft  or  the  ground. 

The  TRAPS  program  was  designed  to  employ  a  combination  of  a 
built-in  1976  standard  atmosphere?  pressure,  temperature,  and  dew¬ 
point  data  profiles  obtained  from  rawinsondes  or  rocketsondes  and 
a  wind  profile  obtained  from  similar  sources.  The  B00MAP2  code 
differs  only  in  that  the  ability  to  process  wind  data  or  nonstan¬ 
dard  atmosphere  has  not  been  verified.  The  meteorological  parame¬ 
ters  play  three  major  roles  in  th<  models  adopted:  a)  specifica¬ 
tion  of  the  atmospheric  pressure  at  the  source  and  key  altitudes, 
b)  defining  the  "effective  speed  of  sound"  iu  the  direction  of 
propagation  as  a  function  of  altitude,  and  c)  defining  the 


effective  speed  of  sound  gradient  seen  by  the  advancing  shc.:k 
wave.  The  consequences  of  these  are  to  displace  the  sonic  boom 
footprint?  affect  the  amplitude  of  the  initial  disturbance  and  to 
determine  the  extent  of  attenuation/enhancement  of  the  shock  wave 
by  the  atmosphere. 

3 . 4  Propagation 

The  sonic  boom  propagation  is  based  upon  the  theory  of 
geometric  acoustics  with  selected  modifications  to  address  its 
peculiar  characteristics.  The  theory  of  geometric  acoustics  is 
valid  when  the  wave  length  is  small  compared  with  characteristic 
macroscopic  scales  of  the  problem.  Such  scales  include  the  radii 
of  curvature  of  the  wave  fronts  and  the  scale  heights  of  the 
atmosphere.  Geometric  acoustics  is  invalid  near  the  aircraft, 
near  a  focus?  and  near  the  boundary  of  and  within  a  shadow  zone. 

In  these  anas  alternative  models  are  required. 

Standard  acoustic  theories  are  linear.  For  sonic  boom 
propagation,  the  cumulative  effect  of  non-linear  effects  over 
large  distance.^  are  significant.  The  cumulative  non-linear 
effects  distort  the  signal  and  produce  shock  waves. 

The  basic  concept  of  the  geometric  theory  is  the  propagation 
of  the  sonic  boom  along  rays,  trajectories  of  points  on  the  wave 
front.  Because  the  wavelength  is  substantially  smaller  than  the 
characteristic  macroscopic  scales  of  the  problem,  it  suffices  to 
trace  only  the  rays  originating  from  the  leading  edge  of  the  shock 
wave.  In  addition,  the  analysis  is  based  on  the  assumption  that 
the  cumulative  non-linear  effects  do  not  affect  the  ray  geometry. 
This  is  an  accepted  assumption  for  most  sonic  boom  problems  of 
interest  including  characterization  of  a  focal  region. 

An  additionally  important  concept  for  the  analysis  is  that  of 
ray  tube.  A  ray  tube  may  be  visualized  as  a  collection  of  rays 
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emitted  from  the  aircraft  initially  displaced  from  each  other  by 
small  times  or  distances.  It  Is  useful  to  define  a  quantity  which 
provides  a  measure  of  the  energy  density  since  it  is  related  to 
the  amplitude  of  the  signal.  Such  a  quantity  is  the  ray  tube  area 
(to  be  defined  more  precisely  later). 

As  a  consequence  of  the  foregoing  assumptions  and  the  hori¬ 
zontally  stratified  atmosphere,  ray  tracing  may  be  performed  using 
a  form  of  Snell's  law.  Using  this  approach,  ray-tube  areas  may  be 
calculated  employing  a  straightforward  numerical  integration.  The 
amplitude  of  the  signal  and  the  amount  of  signal  distortion  may 
then  be  derived  as  a  function  of  these  quantities. 

3.5  Reflection  at  the  Ground 

The  measured  magnitude  of  a  weak  shock  wave  normally  incident 
to  a  perfect  reflector  will  be  twice  the  free  field  overpressures. 
As  a  consequence  of  energy  absorption  by  the  ground,  the  observed 
reflection  factors  are  typically  slightly  less  than  two.  As  the 
angle  of  the  incidence  to  the  ground  becomes  more  oblique,  the 
reflection  factor  decreases,  approaching  one  at  cutoff  (the  loca¬ 
tion  at  which  raypaths  have  been  refracted  to  the  horizontal 
direction  at  the  earth's  surface). 

Although  geometric  acoustics  would  predict  no  sonic  boom  in 
the  shadow  zone  beyond  cutoff,  as  a  consequence  of  diffraction 
effects,  a  low  frequency  rumble  will  be  heard  in  this  zone.  A 
third  phenomenon  that  occurs  near  the  cutoff  boundary  is  the 
focusing  of  the  rays  which  are  turning  up.  As  a  consequence, 
observed  reflectioi  factors  have  the  greatest  spread  near  cutoff. 

The  BQQMAP2  program  models  the  reflection  factor  of  2.0  for 
the  entire  sonic  boom  carpet  out  to  80  percent  of  the  cutoff 
distance.  At  cutoff,  levels  are  reduced  by  10  dB  over  those 
relative  with  a  reflection  factor  of  2.  Beyond  cutoff,  the 
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levels  are  assumed  to  decrease  at  a  rate  of  25  dB  per  decade 
distance  until  a  level  of  80  dB  CSEL  is  reached.  Calculations  for 
a  ray  are  then  terminated.  This  procedure  allows  rapid  truncation 
of  the  boom  levels  beyond  cutoff/  but  avoids  major  discontinuities 
in  noise  level  changes  with  distance. 

3.6  Ray  Tracing 

The  program  assumes  the  atmosphere  (pressures,  temperature, 
and  winds)  to  be  stratified  in  the  vertical  direction,  but  uniform 
in  the  horizontal  direction  and  steady  in  time.  These  assumptions 
impose  stringent  conditions  on  the  possible  paths  of  motion  (rays) 
of  the  wave.  This  motion  is  governed  by  a  variant  of  Snell's  law, 
which  by  virtue  of  the  stratification  of  the  atmosphere,  requires 
the  horizontal  components  of  wave  number,  the  frequency,  and  hence 
the  horizontal  velocity  of  the  phase  surfaces  of  the  wave  to  be 
constant  with  respect  to  the  ground.  This  constant  differs  from 
one  ray  to  another.  When  combined  with  the  requirement  that  the 
net  speed  be  that  of  sound  relative  to  the  air,  it  determines  the 
size  of  the  vertical  component  of  motion,  and  thus,  the  motion 
itself.  The  result  is  that,  for  each  ray,  there  are  combinations 
of  wind  velocity  and  temperature  at  which  it  cannot  exist.  Where 
the  ray  can  exist,  its  path  curves  toward  regions  more  favorable 
to  it;  i.e.,  toward  levels  where  the  sound  speed  is  lower  and/or 
where  the  wind  component  in  its  direction  is  greater  (Figure  5). 
For  each  ray  there  is  a  critical  combination  of  temperature  and 
wind  velocity  that  will  cause  its  vertical  motion  to  slow,  stop, 
and  reverse  (acoustical  cutoff). 


It  should  be  noted  that  a  downward  moving  ray  which  meets  such 
a  reversal  layer  and  turns  away  from  the  ground  will  never, 
because  of  the  stratification  assumptions,  reach  the  ground  no 
matter  what  path  it  subsequently  follows,  but  will  always  reverse 
again  at  the  same  height. 
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FICURE5,.  RAY  CURVATURE  IN  A  STANDARD  ATMOSPHERE 
WITH  NO  IN  DS— HEAD-ON  VIEW. 
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In  a  reference  frame  at  rest  in  the  air  at  the  altitude  of  the 
aircraft  (airborne  reference  frame),  the  normals  to  the  phase  sur¬ 
faces  of  the  wave  can  be  taken  to  have  vector  components  (p,  q,  r) 
in  the  X-,  Y-,  and  vertical  (Z-)  directions,  respectively.  These 
components  represent  the  wave  numbers  in  their  respective  direc¬ 
tions;  the  magnitude  of  this  vector  times  the  sound  speed  is  the 
frequency  (scaled  by  the  aircraft  length),  which  in  the  airborne 
reference  system  is  taken  as  equal  to  the  airspeed  of  the 
aircraft. 

The  tips  of  these  vectors  in  the  airborne  system  must  lie  on  a 
sphere  whose  radius  is  the  aircraft  Mach  number.  In  addition,  it 
can  be  shown  that  the  component  of  the  vector  in  the  direction  of 
the  aircraft  trajectory  must  be  unity.  This  means  that  the  tips 
of  the  vectors  must  lie  in  the  intersection  of  the  Mach-number 
radius  sphere  with  a  plane  normal  to  the  aircraft  motion;  i.e.,  on 
a  circle  which  is  called  the  Mach  circle  (Figure  6).  The  cone 
(ray  cone)  formed  by  the  vectors  from  the  origin  to  the  Mach 
circle  represents  all  the  possible  ray  directions  (in  the  airborne 
reference  system)  from  the  aircraft  at  any  instant;  it3  apex  half- 
angle,  whose  cosine  is  the  inverse  of  the  Mach  number,  is  the 
co-Mach  angle.  An  individual  ray  in  the  cone  i3  specified  by  an 
angle,  ,  which  is  measured  along  the  Mach  circle  from  the 
lowermost  ray,  clockwise  as  seen  by  the  aircraft  pilot  (Figure 
2). 


In  transferring  from  the  airborne  reference  frame  to  one  fixed 
at  the  ground,  the  wave  numbers  p,  q,  and  r  do  not  change.  The 
frequency  “  ,  changes  according  to  the  rule 

+1  vq 

where  u  and  v  denote  the  components  of  the  velocity  difference 
between  the  two  frames  (i.e.,  the  wind  components  at  aircraft 
altitude ) . 


) 
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Because  of  the  stratification  assumptions,  it  may  be  shown 
that  in  any  unaccelerated  reference  frame,  the  parameters  ,  p, 
and  q  do  not  change  as  the  wave  propagates  along  a  ray.  This  is 
the  acoustic  version  of  Snell's  law.  In  addition,  the  following 
relation,  known  as  the  Eiconal  aquation,  holds: 

c2(p2  +  q2  +  r2)  3  {  w  +  up  +  vq ) 2 

at  any  altitude,  where  u,  v,  and  c  are  the  wind  components  and  the 
speed  of  sound,  respectively,  at  that  altitude.  This  defines  an 
admittance  region  in  the  form  of  an  ellipse  for  initial  conditions 
which  can  reach  the  ground.  In  general,  the  admittance  ellipse 
intersects  the  projection  of  the  Mach  circle  at  four  points, 
splitting  the  Mach  circle  into  four  arcs,  two  lying  inside  the 
admittance  ellipse  and  two  outside.  The  arcs  within  the  admit¬ 
tance  region  correspond  to  rays  which  can  penetrate  to  the  ground; 
those  without  cannot.  One  of  the  two  admitted  arcs  will  consist 
predominantly,  or  exclusively,  of  rays  in  the  upper  part  of  the 
Mach  circle  (i.e.,  rising  rays),  the  other  of  the  rays  in  the 
lower  part  (descending  rays). 

The  program  contains  a  routine  to  carry  out  the  construction 
of  the  admittance  ellipse  at  ground  level,  and  to  determine  for 
which  initial  orientations,  <j>  ,  the  rays  will  lie  inside  the 
ellipse.  Since  interest  is  in  the  sonic  boom  on  the  ground,  the 
program  declines  to  trace  rays  outside  the  admittance  ellipse 
which  saves  computer  time  and  printout. 

For  a  ray  within  the  admittance  ellipse  initial  conditions  for 
ray  tracing  are  determined  a3  follows.  The  initial  position  of 
the  ray  is  taken  at  the  tip  of  the  aircraft  (x(t^),  yft^), 
z(t&))  with  the  initial  direction  of  propagation  given  by  the 
ray  bank  angle,  <f>  ,  the  aircraft  Mach  number,  M,  the  aircraft 
heading,  e  ,  and  climb  angle,  y  (both  relative  to  the  wind).  The 
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initial  wave  mmbers  and  frequency  are  then  calculated  as 


2  2 

Recurvature  (ray  direction  reversal)  occurs  when  (A/c)  <  p 

2  ® 

+  .  In  order  to  forecast  recurvature  and  the  need  for  smaller 

step  size  the  following  expression  is  used  for  the  z  ray 

"acceleration" . 

[p0  u,  +  q0  vz 

3.7  :%ay  Tube  Area 


-  c  c. 


t 


+q 


At  each  instant  of  supersonic  flight,  the  aircraft  emits  a 
cone  of  rays,  each  of  which  is  singled  out  by  specifying  an  angle, 
<t»  .  The  set  of  rays  which  leave  the  aircraft  at  neighboring  times, 
between  t^  and  t&  +  At^,  and  at  neighboring  angles,  between 
(jj  and  $  +  A<p  ,  form  a  ray  tube. 

[  The  total  acoustic  energy  in  a  ray  tube  has  been  shown  to  be 

constant  (for  linear,  inviscid  processes)  by  Blokhintzev  who 
i  formulated  an  invariant  relating  the  ray  tube  area  to  the  inverse 

square  of  the  amplitude.  This  invariant  is  used  in  both  the  TRAPS 
and  the  B00MAP2  programs. 


The  programs  define  the  ray  tube  area  as  neither  a  horizontal 
section  nor  a  cross  section,  but  as  a  section  cut  by  the  wave 
phase  surfaces  within  a  unit  time,  i.e„,  a  section  normal  to  the 
wave  normals. 


Defined  in  this  way,  the  ray  tube  area  is  always  finite,  and 
is  a  Galilean  invariant.  That  is,  it  is  a  quantity  whose  value 
does  not  change  when  measured  by  an  observer  moving  at  any 
constant  velocity.  Since  the  amplitude  of  the  sonic  boom  is 
clearly  a  Galilean  invariant,  as  are  the  pressures,  temperatures, 
densities,  sound  speeds,  and  other  physical  quantities  in  the 
Blokhintzev  invariant,  this  definition  is  the  most  appropriate. 
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The  program  computes  the  area  as  a  determinant  (called  a 
Jacobian)  formed  from  partial  derivatives  of  coordinates  with 
respect  to  the  ray  parameters  $  and  t^.  These  partial  deriva¬ 
tives  are  found  by  integrating  equations  similar  to  the  ones  used 
to  track  the  rays,  and  are  in  fact  derived  from  them  as  will  be 
discussed  subsequently. 

This  technique  is  better  than  the  alternative  of  actually 
tracking  neighboring  rays  and  computing  the  area  of  the  figure 
formed  by  the  endpoints,  since  over  the  distances,  even  rays  which 
are  initially  very  close  can  spread  over  considerable  distances. 
Furthermore,  area  computations  of  that  type  are  so  sensitive  to 
round-off  errors  in  position  that  the  error  may  be  many  times  the 
actual  area. 

Specifically,  the  ray  tube  area  is  defined  as  the  Jacobian,  J, 
given  by 

-  _  9  (x  ,y ,  z) 

J  “  3(Y,t a,4>) 

where  4<  is  the  phase 

(Ray  tube  area  has  units  of  length  squared  per  unit  time.)  The 
Blokhintzev  invariant  alluded  to  earlier  is  3imply 

J  p  2 
*o 

- *  constant. 

pa  6  c 

V 

Along  any  ray  the  Jacobian  value  may  be  calculated  by  taking  the 
partial  derivatives  of  the  ray  tracing  equations  with  respect  to 
4'/  ta,  and  $  and  interchanging  the  order  of  differentiation. 

Since  the  ray  bank  angle,  4> ,  does  not  influence  the  position 
of  the  aircraft,  ax/s^,  3 t/a 45  and  32/341  are  initially  zero.  The 
remaining  initial  conditions  are  found  by  differentiating  the  wave 
number  initial  conditions  with  respect  to  $ 
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Similarly,  initial  conditions  are  obtained  on  the  derivatives 
of  x,  y,  z.  A,  to,  p,  q,  and  r  with  respect  to  t^.  To  differen¬ 
tiate  with  respect  to  t^  with  Y  fixed,  the  value  of  dx/dt  along 
the  ray,  i.e.  u  +  cp  must  be  subtracted  from  x'tt^).  Similar 
adjustments  must  be  made  for  the  other  coordinates,  y  through  r. 

Since  the  third  parameter  defining  tne  rays  is  the  phase  Y 
(which  is  identified  with  position  5  behind  the  nose  and  has  the 
unusual  dimension  of  length),  initial  conditions  on  derivatives 
with  respect  to  it  are  not  required.  Rather,  since  the  vectors 
3  (x,y,z  )/3<t>  and  3  (x  y,z)/3tA  are  both  expected  to  lie  in  the 
surface  Y*  const.,  the  vector  3  (x,y  ,z  )/‘34'  is  replaced  with  the 
normal  vector  VY  /  j  VYp  or  -  (p , q,  r  )/ (p2+q2+r2  )  without 
changing  the  Jacobian. 

The  Jacobian  technique  leads  to  a  ray  tube  area  that  varies  in 
a  continuous  manner  as  the  ray  is  traced,  and  even  the  rate  of 
change  of  area  with  position  along  the  ray  is  continuous  so  long 
as  the  gradients  of  wind  and  sound  speed  are  continuous  in  the 
atmosphere  model.  Where  the  gradients  are  discontinuous  (and  this 
occurs  at  each  height  at  which  either  tempera!  ure  or  wind  is  input 
or  taken  from  the  Standard  Atmosphere),  the  rate  of  change  of  area 
(but  not  the  area)  undergoes  a  jump.  The  amount  of  this  jump  is  a 
continuous  function  of  the  ray  normals,  which  are  themselves 
continous. 
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The  result  of  this  is  that  the  ray  tube  areas  on  the  ground 
and  the  amplitudes  are  continuous  functions  of  the  ray  parameters. 
Except  when  the  ray  tube  area  is  zero,  or  at  the  edge  of  the 
carpet,  they  are  also  continuous  functions  of  position  on  the 
ground. 

3.8  Signal  Propagation 
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In  the  linearized  acoustic  theory,  the  wave  form  of  the 
pressure  travels  along  the  ray  unchanged  except  for  amplitude 
changes  governed  by  the  Blokhintzev  invariant.  At  least  below  the 
mesopause,  effects  of  viscosity  and  heat  conduction  are  too  small 
to  seriously  affect  this  concept. 

Pressure  waves  of  this  amplitude  are  governed  by  a  non-linear 
theory,  and  although  the  non-linear  effects  are  small  over  any 
given  region  up  to  seme  tens  of  wavelengths  in  size,  they  do 
accumulate  and  are  responsible  for  the  typical  N-wave  profile  of 
the  direct  sonic  booms  and  the  bulk  of  dissipitation  of  acoustic 
energy  between  the  aircraft  and  ground. 

In  terms  of  supersonic  flow,  the  sonic  boom  is  "weak,"  and  the 
program  applies  a  weak  shock  tube  theory  due  to  G.B.  Whitham  to 
the  propagation  of  the  3onic  boom  in  ray  tubes.  In  general,  an 
overpressure  at  a  given  point  in  the  wave  form  so  increases  the 
air  speed  and  sound  speed  at  its  location  that  it  seems  to  over¬ 
take  a  lesser  overpressure  located  ahead  of  it.  The  amount  of  the 
overtaking  is  governed  by  a  quantity  termed  the  age,  which 
increases  along  a  ray  at  a  rate  proportional  to  the  amplitude,  and 
inversely  proportional,  among  other  terms,  to  the  square  root  of 
the  ambient  air  density.  The  age  is  given  by  the  expression 

£ 

Alt)  =  .5  (T+l)  ]  y/2,^ 

t  ,2  ,,*5 

a  (c  PaJ) 
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(Age  has  units  of  length  multiplied  by  time  divided  by  the  square 
root  of  mass. ) 

The  shift  in  phase  may  be  expressed  as 


H"  =  V  +  Ap  i  ( <f> ,  4* ) 
where 

Pi(<t.H')  =  cA  /p^  M  F  (<J> ,  H') 

(F  has  units  of  length*5  ,  pj_  has  units  of  mass*5  per  unit  time.) 

When  a  section  of  the  waveform  actually  overtakes  one  ahead  of 
it,  the  choice  among  the  three  or  more  possible  values  of  over¬ 
pressures  is  resolved  by  fitting  a  shock  (pressure  jump),  thereby 
cutting  off  the  lobes  of  the  overtaking  and  overtaken  portions. 

To  conserve  mass,  the  shocks  are  so  placed  as  to  balance  the  area 
within  the  cutoff  lobes  using  the  so-called  "equal  area  rule" 

( Figure  7  ) . 

The  pressure  far  field  signature  is  calculated  in  terns  of  the 
above  quantities  as 

p(<j),4')  =  c  (Apa)*s  j  15  pi($,'n 

When  the  ray  tube  area  reverses  sign  along  a  ray  the  geometric 
theory  implies  that  the  pressure  is  infinite  at  the  point  where 
the  area  is  zero.  Thi3  point  is  called  a  caustic  point.  Infinite 
pressures  are,  of  course,  contrary  to  reality  and  are  a  conse¬ 
quence  of  the  failure  of  this  theory. 

In  fields  other  than  acoustics,  such  as  water  wave  theory  or 
optics  to  which  ray  theory  applies,  a  more  general  theory  known 
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FIGURE  7  .  SIGNATURE  AGING  PROCESS  ADAPTED 
FkOM  HAYtS  ET  AL  (1969), 


a3  Uniform  Asymptotic  theory  may  be  applied.  Indeed,  this  more 
general  theory  holds  for  linearized  acoustics  as  well,  and  can  be 
used  to  determine  the  shape  of  the  wave  departing  the  caustic, 
given  the  shape  of  the  wave  approaching  the  caustic.  After  pass¬ 
age  through  the  caustic,  ordinary  ray  theory  holds  once  again  and 
the  program  may  resume,  now  propagating  the  new  signature.  It  is 
t'  s  technique  which  the  B0QMAF2  program  uses  to  continue  the 
e .  jlution  of  the  sonic  boom. 

It  is  a  conclusion  of  the  Uniform  Asymptotic  theory,  to 
whichever  physical  procesr.  it  has  been  applied,  that  the  Fourier 
components  of  the  outgoing  signal  are  the  same  as  they  would  be 
expected  to  be  from  the  naive  ray  theory,  except  that  each  one  has 
been  shifted  forward  one  quarter  wavelength.  Since  the  shorter 
wavelength  components  advance  less  than  the  longer  components,  the 
shape  of  any  complex  waveform  can  change  significantly. 

This  transformation  is  commonly  known  by  the  name  of  "90 
degree  phase  shift"  (since  there  are  360  degrees  in  a  full  wave 
cycle).  Hence,  there  is  a  temptation  to  perform  it  by  actually 
taking  a  finite  Fourier  transform,  changing  the  coefficients,  and 
inverting.  However,  even  with  the  Fast  Fourier  Transform,  this  is 
an  extremely  inefficient  procedure. 

The  reason  lies  in  the  shape  of  the  input  signal,  which  by  the 
time  of  caustic  passage  has  usually  aged  into  a  nearly  N-wave 
form.  The  transform  of  the  N-wave  has  two  very  thin  peaks •{ logar¬ 
ithmic  discontinuities)  located  where  the  jumps  were  (Figure  8). 

To  resolve  these  peaks  requires  a  number  of  very  closely  spaced 
points  in  their  immediate  vicinity.  Elsewhere,  the  waveforms  are 
smooth  and  such  close  spacing  in  extremely  wasteful  of  computer 
resource.  In  particular,  a  much  wider  spacing  should  be  used  far 
ahead  of  and  far  behind  the  original  waveform.  But  finite  Fourier 
transforms  require  a  uniform  spacing  of  points,  forcing  a  choice 
between  inadequate  resolution  and  waste  of  resources. 
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The  program  uses  an  alternative  to  the  above  Fourier 
techniques,  called  the  Hilbert  Transform.  This  is  an  integral 
transform  with  a  singular  kernel  whose  Fourier  equivalent  happens 
to  be  the  90  degree  phase  shift;  it  has  the  advantage  that  it  may 
be  evaluated  at  an  arbitrary  selection  of  points  whose  spacing  may 
be  chosen  with  the  above  principles  in  mind. 

In  the  program,  the  sonic  boom  signature  is  taken  through  the 
following  evolutionary  steps; 

(i)  Compute  the  age  until  the  ground  or  a  caustic  is 
encountered  ? 

(ii)  Age  the  signature  and  fit  shocks  as  appropriate; 

(iii)  If  at  a  caustic,  perform  the  Hilbert  Transform  and  create 
a  new  signature; 

(iv)  Continue  with  step  (i)  until  final  ground  contact. 

In  this  analysis,  the  Uniform  Asymptotic  theory  must  be 
regarded  as  an  approximation  in  that  the  shocks  of  the  N-wave 
indicate  the  operation  of  non-linear  effects,  and  the  theory 
applies  to  linear  systems  which  is  reinforced  by  the  appearance  of 
infinities  in  the  Hilbert  Transform  of  the  N-wave.  In  reality, 
however,  the  N-wave  with  the  shocks  is  an  approximation  to  the 
actual  signature.  Since  the  sonic  boom  is  weak,  in  the  3ense  of 
supersonic  flow  theory,  the  shocks  are  not  strong,  well 
established  features.  Measurements  often  show  a  "rise  time"  for 
the  shocks  of  between  1/30  and  1/10  of  the  length  of  the  N-wave, 
presumably  due  to  3ome  form  of  turbulent  dispersion.  With  such  a 
"thick  shock,"  the  infinities  in  the  wave  form  all  disappear,  and 
the  Uniform  Asymptotic  theory,  if  carried  out,  would  lead  to 
finite  overpressures  up  to  and  past  the  caustic  surface.  This 
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result  places  the  validation  of  the  Uniform  Asymptotic  theory  on 
the  same  level  as  ray  theory,  as  an  approximation  to  the  linear 
acoustic  equations,  and  the  validation  of  the  linear  theory  as  an 
approximation  to  the  non-linear  theory  on  the  same  level  near  the 
caustic  as  elsewhere. 

If  a  caustic  point  lies  near  the  ground,  it  is  important  to 
characterize  the  sonic  boom  signature  at  this  location.  There  are 
three  types  of  caustic  points:  a  smooth  caustic,  a  cusped  caus¬ 
tic,  and  a  perfect  focus.  The  smooth  caustic  lies  along  a  surface 
containing  continuous  focusing  for  a  range  of  initial  times  and 
ray  angles.  The  ground  intersection  is  a  line.  A  cusped  caustic 
has  an  infinitesimal  perfect  focus  along  a  curve  intersecting  the 
ground  in  a  point.  Similarly  a  perfect  focus  will  (at  most) 
intersect  the  ground  in  a  point;  it  esults  from  a  finite  wave 
element  focusing  to  a  point.  The  program  includes  a  model  of  the 
smooth  caustic.  It  is  the  most  frequent  type  and  affects 
significantly  larger  areas  than  the  other  two  types. 

The  signature  calculated  for  a  smooth  caustic  is  based  on  the 
smooth  caustic  similitude  solution  developed  by  Seebass,  This 
solution  requires  determining  the  relative  curvature  of  the  caus¬ 
tic  surface  in  the  direction  of  collapse  with  respect  to  the  ray 
curvature.  Th  calculations  proceed  as  follows: 

Auxiliary  rays  are  traced  to  determine  the  direction  of  ray 
tube  collapse.  The  ray  along  which  focusing  was  detected  may  be 
characterized  by  the  time,  t0/  and  the  ray  bank  angle,  <f>  o,  at 
which  it  was  emitted.  The  auxiliary  rays  are  perturbed  from  the 
reference  ray  as  follows: 
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where 

A<t  -  sine  (<|>)x  0.5  degrees  *  sine  ($)  x  8.73~3  x  10  radians 

4.6  ACL 
A<j>  c 

At  »  - - -  seconds 


where  ACL  =  aircraft  length 

-  speed  of  sound  at  the  aircraft 

At  the  focus  the  vectors  from  ray  1  to  ray  3  and  ray  3  to  ray 
4  are  computed  and  designated  t  13  and  j 34  respectively. 

Additional  focusing  ray  tubes  are  calculated  as  follows: 
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Summary  of  Raytubes  to  be  Traced 


i 


Origin 


Focus 


Tube 

4> 

t 

Location 

Time 

1. 

Original  Focusing 

Tube 

<P 

-► 

po 

T1 

2. 

First  Auxiliary 

♦l 

fcl 

m  fco 

AS 

-  mca 

"P0  "  As 

feet" 

T2 

3. 

Second  Auxiliary 

*2 

t2 

*  fco 

+  ££- 
Mca 

"P0  +  As 

feet" 

T3 

4. 

Third  Auxiliary 

$3 

«  t 

o 

2  As 

‘  Mca 

"Po  -  2As 

feet" 

T4 

As  *  800  feet 

- 

‘j  *  ‘o  +  *inc  vj  (tj  '  ‘o’ 

v .  *  average  aircraft  velocity  between  t  and  t . 


ip . 

xnc 


^13  *  ^34  A 4> 

71 - 2~  vAt 

I  A34 


v  =  aircraft  speed  at  time  t 


The  focus  times  for  the  four  focusing  ray  tube  are  examined  to 
assure  that  a  smooth  focus  is  being  treated:  T4<T2<Tl<T3 '  n0t 
the  focus  is  discarded  (not  a  smooth  caustic). 


The  points  along  the  reference  ray  upstream  of  the  focus  are 
fitted  with  a  circular  arc  to  estimate  the  ray  curvature  at  the 
focus,  kr.  Similarly,  the  caustic  curvature  <c;  is  calculated  from 
the  first,  three  focusing  ray  tubes. 
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The  relative  curvature  vector  upon  which  the  amplitude  of  the 
focus  overpressure  is  calculated  as 
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and  the  relative  radius  of  curvature  is  given  by 
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R 


rel 


•rel 


In  order  to  proceed  with  the  calculations,  the  program  now 
searches  for  a  point  on  the  reference  ray  for  which  the  peak 
overpressure  matches  the  peak  calculated  from  the  focus  solution. 
This  is  implemented  in  an  iterative  process  as  follows.  At  a 
point  sufficiently  upstream  of  the  focus  the  peak  overpressure  is 
calculated.  The  pressure  coefficient,  Cp,  for  the  largest 
overpressure  is  calculated  as 


'P  = 


A  P 


2  r  p 


The  distance  from  the  caustic  along  the  ray  to  this  reference 
point,  S,  is  calculated  by  summing  the  arc  segments.  The  distance 
normal  to  the  focus  is  estimated  as 


.  52rel 

This  distance  is  now  used  to  estimate  focal  zone  boundaries  (where 
the  incoming  signal  matches  the  peak  focal  overpressure)  first  in 
the  direction  normal  to  the  focuss 
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and  then  in  a  tangential  direction 


X*  -  ^  RrelY‘ 


If  the  selected  point  is  more  distant  than  X*,  this  procedure 
is  repeated  until  a  point  at  this  distance  is  located  (by  interpo¬ 
lation  as  necessary).  A  detailed  signature  is  then  developed 
using  the  full  Gill/Seebass  solution  for  each  shock  in  the 
solution.  If  the  ground  level  lies  between  the  focus  and  a  dis¬ 
tance  X*  upstream  of  the  focus,  the  focus  signature  (scaled  by  the 
appropriate  reflection  factor)  is  taken  as  the  ground  level 
signature. 


If  the  ground  level  is  sufficiently  far  downstream  of  the 
f o»_i s  the  TRAPS  postfocus  solution  discussed  earlier  provides  an 
adequate  representation  of  the  ground  overpressures.  In  order  to 
assure  that  the  ground  level  is  sufficiently  far  downstream  of  the 
focus,  the  same  criteria  are  employed  as  were  used  on  the  upstream 
side  of  the  focus.  The  pressure  coefficient  corresponding  to  the 
peak  pressure  in  the  TRAPS  focu3  solution  i3  compared  with  the 
pressure  coefficient  from  the  peak  focus  overpressure.  If  the 
focus  solution  is  larger  the  TRAPS  postfocus  solution  is  accepted. 
Otherwise,  the  focus  solution  is  used  as  the  ground  level  free 
field  signature. 


3 . 9  Implementation  of  Ray  Tracing  In  B0QMAP2 


At  any  given  aircraft  track  time  t,  the  aircraft  position  and 
velocity  vectors  are  known,  and  the  acceleration  vector  and  inter- 
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at  ground  level,  defined  in  Section  3.6,  determines  the  range  of 
the  initial  ray  bank  anple  which  will  reach  the  ground.  Starting 
with  the  largest  negative  value,  rays  are  traced  at  increments 
of  1°  if  the  aircraft  is  above  15,000  ft  and  2°  if 
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the  aircraft  is  below  15,000  ft,  until  the  largest  positive  value 
of  $  is  reached. 


If  during  the  tracing  of  a  ray,  a  caustic  is  encountered  in 
the  "ground  zone",  defined  as  the  region  1000  ft  above  the  ground 
to  1500  ft  below  the  ground,  the  ray  history  and  caustic  location 
are  stored,  the  ray  tracing  terminated  and  the  next  value  of  <j>  is 
traced.  After  all  possible  Rvalues  have  been  traced  for  the  time 
t,  the  stored  caustic  locations  are  used  to  estimate  no  more  than 
two  ray  bank  angles4>,  at  which  the  caustic  surface  either  crosses 
the  ground  or  is  closest  to  it. 

Startinq  at  these  estimated  $  values,  rays  are  traced  at  $ 
intervals  of  +0.1°  or  +0.2°  (depending  on  aircraft  altitude)  until 
the  caustic  is  no  longer  in  the  "ground  zone".  Using  the  theory 
outlined  in  Section  3.8  the  focal  zone  width  is  estimated  for  an 
angle  $  ,  If  the  ground  is  within  the  focal  zone,  the  focus 
overpressure  and  signatures  are  used,  otherwise  the  original  TRAPS 
solution  is  used. 


The  relative  curvature  of  the  ray  and  the  caustic  surface  are 
used  in  calculating  the  focal  width.  The  tracing  of  the  auxiliary 
ray  tubes  to  define  the  caustic  curvature  necessitate  interpola¬ 
tion  beween  at  least  three  adjacent  aircraft  track  times.  If  the 
original  aircraft  position  and  velocity  vectors  are  "noisy",  the 
acceleration  and  jerk  vectors  (even  after  smoothing)  can  show 
major  variations  with  time,  and  hence,  the  caustic  surface  can  be 
irregular  or  form  cusps.  The  caustic  curvature  may  show  major 
variations  with  <j>  angle  and  give  unrealistic  large  values  of  focal 
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zone  width  extended  more  than  2500  ft  above  the  ground,  the  ray 
was  rejected.  In  addition,  the  overpressures  calculated  on  the 
ground  for  all  <j>  values  (both  focus  and  TRAPS  solutions)  are 
checked  for  singularities  which  give  too  large  a  variation  with,),, 
and  these  rays  are  rejected. 
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3.10  Calculation  of  CSEL  from  Signatures 


Typically,  there  is  a  simple  relationship  between  over¬ 
pressure  and  C-weighted  Sound  Exposure  Level  (CSEL)  for  sonic  boom 
signatures,  which  are  sufficiently  far  away  from  a  focus  (Ref.  9) 

lCE  3  kPK  -  26 

However,  in  the  focal  and  post-focal  zones,  the  pressure  signature 
changes  radically  (Figure  9)  and  a  more  accurate  estimate  of  CSEL 
may  be  required,  which  is  done  by  performing  a  Fast  Fourier 
Transform  (FFT)  of  the  signature,  applying  a  C-weighted  filter  to 
the  spectrum  and  integrating  to  give  CSEL. 

The  signature,  produced  by  B00MAP2 ,  gives  pressures  (including 
shocks)  at  irregularly  spaced  time  intervals.  Based  upon  the 
parametric  study  of  waveforms  in  Appendix  C,  Ref.  9, 

(a)  the  effect  of  the  rise  time  of  the  shocks  on  CSEL  is 
insignificant; 

(b)  A  Nyquist  frequency  of  1000  Hz  is  adequate  (i.e.  time 
interval  =0.5  milliseconds). 


The  BOOMAP2  signature  is  therefore  modified  by  first 
identifying  the  shocks  and  separating  them  by  0.5  milliseconds. 

The  time  of  maximum  overpressure  is  used  to  define  a  new  time  var¬ 
iable,  spaced  at  equal  intervals  of  0.5  milliseconds,  and  the 
pressure  signature  is  interpolated  at  these  times.  An  FFT  is  per¬ 
formed,  after  extending  the  signature  with  zeros  to  the  necessary 
(power  of  2)  number  of  points.  An  analytical  C-weighted  filter  is 
applied  to  the  spectrum  and  the  spectrum  is  then  integrated  to 
give  the  CSEL  value. 


3.11  Calculation  of  Scratchpad  Contours 


For  the  scratchpad  contour  plots  to  be  generated,  the  user 
must  first  select  the  3orties  to  be  processed.  The  BOOMAP2  pro¬ 
gram  then  accesses  the  database  containing  the  rays  associated 
with  the  selected  sorties.  All  rays  emitted  from  the  aircraft 
at  the  same  time  are  then  read  from  the  database  and  sorted  by 
angle  This  processing  continues  until  all  the  rays  for  a 
selected  sortie  are  processed,  or  a  4-1/2  second  time  gap  is 
found.  This  is  known  as  a  flight  segment. 

If  no  caustic  rays  are  found  in  the  flight  segment,  then  a 
scratchpad  contour  plot  will  not  be  generated  for  that  flight 
segment.  The  flight  segment  is  then  appended  to  a  temporary 
file  only  if  it  contains  caustic  rays.  This  process  is  repeated 
until  all  the  selected  sorties  are  processed. 

After  the  selected  sorties  are  processed  and  the  temporary 
file  is  created,  the  B00MAP2  program  then  reads  an  entire  flight 
segment.  It  then  converts  the  pressure  in  pascals  to  either  the 
maximum  overpressure  in  psf  or  the  maximum  sound  pressure  level 
(SPL)  depending  upon  which  metric  the  user  specified.  At  the 
same  time  the  maximum  pressure  for  the  entire  flight  segment  is 
found.  Based  on  this,  ten  contour  levels  are  selected  either  in 
dB  or  psf.  If  the  maximum  SPL  has  been  selected,  the  contour 
levels  are  in  intervals  of  2.5  dB.  If  overpressure  in  psf  has 
been  selected,  the  contour  levels  are  in  a  ratio  of  two  to  one 
(0.25,  0.5.  1.  2.  ...). 

The  next  step  divides  the  flight  segment  into  sections 
based  on  initialization  time  from,  the  aircraft.  These  sections 
are  then  searched  for  the  selected  contour  level.  If  the 
selected  contour  level  is  found,  a  linear  interpolation  is  done 
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to  find  the  x,  y  location  of  that  pressure;  up  to  four  x,  y 
points  may  be  returned  for  each  time  slice.  If  no  points  are 
found  for  the  flight  segment,  a  search  for  the  next  contour 
level  begins.  This  processing  continues  until  all  ten  contour 
levels  have  been  searched  or  points  have  been  found  in  at  least 
two  different  time  slices  for  the  flight  segment.  Then,  the 
points  that  have  been  identified  are  sorted  on  termination  time 
in  order  to  allow  the  outer  edges  of  the  contour  to  be  connected 
in  increasing  x  and  y  ground  locations.  Once  the  points  are 
connected,  the  contour  is  plotted  on  the  scratchpad.  If  possi¬ 
ble,  three  contour  levels  are  plotted  on  each  scratchpad.  This 
processing  continues  until  all  selected  sorties  have  been  pro¬ 
cessed. 

3.12  Calculation  of  Average  Overpressure  Values  for  GPCP 

Contouring  Program  Processing 

The  ray  tracing  code  (see  Section  3.9)  results  in  calculat¬ 
ed  overpressure  values  at  an  irregularly-spaced  array  cf  ground 
positions.  The  GPCP  contouring  program  requires  as  input  that 
the  overpressures  be  defined  at  grid  points  that  are  equally 
spaced  in  a  rectangular  pattern.  Thus,  it  is  necessary  to  de¬ 
velop  algorithms  to  calculate  values  at  the  GPCP  grid  points 
that  are  based  on  the  calculated  ray  tracing  overpressures  in 
the  vicinity  of  each  grid  point.  Each  grid  point  may  be  influ¬ 
enced  by  none,  one,  or  many  calculated  ray  values. 

The  gridding  algorithm  for  determining  values  at  the  GPCP 
grid  points  requires  three  arrays;  a  master  array,  an 
accumulator  array  and  a  counter  array.  Each  of  these  arrays  is 
dimensioned  102  by  102  with  each  cell  representing  a  grid  point. 
The  spacing  between  grid  points  is  set  at  2500  feet.  The  combat 
training  range  geographic  center  is  located  in  the  center  of 
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each  array  allowing  125,000  ft  in  each  direction  from  the  range 
center. 

The  first  step  in  the  gridding  algorithm  is  to  initialize 
the  master  grid  to  a  value  of  1.0  pascals.  This  reduces  the 
extent  of  the  steep  contours  resulting  from  an  otherwise  zero 
background . 

The  rays  that  coincide  with  the  selected  sortie  are  then 
read  from  the  database,  generated  by  the  B00MAP2  program,  one 
sortie  at  a  time.  The  data  is  then  sorted  by  $  angle.  If  there 
are  any  gaps  in  the  data  greater  than  two  degrees,  a  linear 
interpolation  is  used  to  fill  the  gap  in  one  degree  increments. 
The  rays  for  the  selected  sortie  are  then  sorted  on  termination 
time . 


These  rays  are  then  processed  in  5.5  second  time  slices. 
Each  ray  in  the  time  slice  is  mapped  to  the  four  closest  grid 
points  to  where  it  terminates.  The  accumulator  array  is  then 
updated  by  adding  the  squared  pressure  of  the  ray,  in  pascals, 
to  the  current  value  of  the  grid  point.  The  four  corresponding 
counter  array  grid  points  are  then  incremented  by  a  value  of 
one.  This  processing  continues  until  all  the  rays  for  the  time 
slice  have  been  processed. 

After  5.5  seconds  have  elapsed,  the  master  grid  array  is 
then  updated  by  dividing  the  accumulator  array  by  the  counter 
array  and  adding  it  to  the  current  value  of  the  master  array. 
The  accumulator  array  and  the  counter  array  are  then  zeroed. 
This  processing  is  continued  until  all  the  rays  for  a  selected 
sortie  are  processed. 
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After  the  selected  sortie  has  been  processed,  the  above 
process  is  repeated  until  all  of  the  user's  selected  sorties 
have  been  processed.  Once  this  is  done,  the  master  grid  is  then 
divided  by  the  number  of  supersonic  sorties  and  the  square  root 
is  taken  to  get  the  KMS  value.  The  result  i.s  converted  to  CSEL 
values  using  the  following  formula,  (20  *  loglO(n)  +  68).  At 
this  point,  the  master  grid  values  can  be  fed  into  the  GPCP  con¬ 
touring  program. 
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4.0  AIRCRAFT  F-FUNCTIONS  FOR  TBS  MODIFIED  TRAPS  PROGRAM 


As  discussed  in  Section  2,  a  sonic  boom  propagation  program  | 

must  be  provided  with  a  description  of  the  disturbance  of  the  I 

atmosphere  generated  by  the  supersonic  aircraft.  It  may  be 
described  in  terms  of  the  overpressure  wave  form  that  an  observer 
near  the  aircraft  would  measure  as  the  aircraft  flies  by  or,  j 

equivalently,  in  terms  of  a  theoretically  derived  F-function 
defining  the  flow  near  the  aircraft.  Both  the  original  TRAPS 
program  and  the  B00MAP2  computer  program  use  the  latter  approach. 

The  F-function  is  developed  from  consideration  of  the  aircraft 
geometry  and  its  lift  distribution  (see  more  detail  in  Section  2). 

The  F-function  can  generally  be  separated  into  area  and  lift 
components.  The  original  TRAPS  program  provided  for  the 
F-function  described  in  terms  of  these  two  components.  However, 
the  BOOMAP 2  code  employs  a  simplified  F-function  based  on  a  method 
developed  by  Carlson  (Reference  4).  This  F-function  can  be  char¬ 
acterized  as  an  N-wave  with  a  rise  time  much  shorter  than  the  dur¬ 
ation  of  the  entire  disturbance.  The  peak  amplitude  of  the  I 

F-function  is  then  taken  as  a  function  of  the  aircraft  shape  ! 

i 

factor,  Ks  shown  below:  | 


F-function 
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In  the  modified  TRAPS  program  implementation,  the  F-function 
is  characterized  as  having  a  rise  time  l/100th  the  duration  of  the 
entire  disturbance.  However,  rise  time  is  not  critical  in  the 
overpressure  calculations  since  it  is  the  area  of  the  F-function 
displayed  that  is  important  and  this  area  is  independent  of  the 
rise  time. 


The  approach  used  to  develop  simplified  F-functions  follows 
that  of  Carlson.  Thus,  Ks  is  estimated  from  curves  showing  Ks 
as  a  function  of  (using  Figure  4  of  Reference  4).  Kl  is 
defined  as  follows: 
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Based  on  thi3  approach,  the  BQ0MAP2  program  has  the 
F-functions  as  shown  in  Table  1.  It  lists  the  aircraft  in  the 
current  program  together  with  typical  aircraft  lengths  and 
weights,  together  with  typical  KL  and  Ks  functions. 


F-functions  for  new  aircraft  can  be  developed  following  the 
procedures  of  Carlson,  but  these  procedures  are  a  simplification 
of  actual  situations,  because  the  F-functions  change  with  aircraft 
accelerations  and  also  vary  with  altitude  and  speed.  However, 
these  changes  generally  are  small  compared  with  other 
uncertainties  in  the  calculations,  and  are  not  critical  unless  one 
is  concerned  with  developing  a  detailed  calculation  for  specific 
test  conditions,  a  situation  beyond  the  routine  application  of  the 
B00MAP2  program. 
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TABLE  1 


TYPICAL  AIRCRAFT  LIFT  PARAMETER  AND 
SHAPE  FACTOR  VALUES  FOR  SONIC  BOOM 
CALCULATIONS 


AIRCRAFT 

LENGTH 

FT 

WEIGHT 

KLBS 

TYPICAL 

KL 

TYPICAL 

KS 

B~  IB 

147.0 

453 . 0 

0.0060 

0 . 0910 

F-4 

50. 2 

56 . 0 

0 . 0040 

0 . 0860 

F-5 

46.6 

19.6 

0.0025 

0 . 0642 

F-0 

54.5 

32.3 

0.0035 

0 .0070 

F- 1 4 

62.7 

56.7 

0.0040 

0 . 0873 

F- 1 5 

63.0 

42 . 3 

0.0030 

0 . 0638 

F-16 

47.6 

23.3 

0.0030 

0 . 0838 

F-18 

56.0 

49.3 

0.0050 

6.0900 

F-20 

46.5 

26.1 

0.0035 

0 . 0643 

F-i  01 

71 . 1 

48.4 

0.0030 

0 . 0860 

F-104 

54.8 

21 .4 

0.0025 

0 . 0642 

F- 105 

64.2 

42.7 

0 . 0030 

0.0069 

F-106 

70.8 

34 . 2 

0.0020 

0 . 0840 

F-lll 

75.5 

95.0 

0.0050 

0 . 0692 

SR-71 

107.4 

161 . 0 

0.0100 

0 . 0370 

T-3B 

46.3 

11.2 

0 . 0020 

0.0642 

5.0  SELECTION  OE  THE  COMPUTATIONAL  PROGRAM  TRAPS  FOR  SONIC  BOOM 
CALCULATIONS 

The  calculation  for  sonic  booms  for  general  maneuvers  in  a 
real  atmosphere  is  sufficiently  complex  to  require  a  computerized 
model.  A  number  of  computer  models  exist,  all  of  which  rest  on 
the  identical  theory  for  nonfocusing  cases.  Major  differences  lie 
in  the  computational  philosophy  and  added  features.*  Four  exist¬ 
ing  models  were  considered  for  the  current  programs 

A.  SABER  (Ref.  12)  developed  by  the  J.  H.  Wiggins  Company  for 
USAF  WSMC ,  which  is  descended  from  the  Thomas  program. 

B.  TRAPS  (Ref.  3)  developed  by  the  NOAA  Air  Research  Laboratory 
has  its  philosphical  and  conceptual  origins  in  the  ARAP 
program. 

C.  FOBOOM  (Ref.  13)  developed  by  Wyle  Laboratories  is  an  exten¬ 
sion  of  the  THOMAS  program  which  can  compute  boom  signatures 
at  focal  zones. 

D.  SABERII  (Ref.  14),  an  evolutionary  development  of  the 
WYLE/MSFC  model  for  applications  to  Space  Shuttle  ascent. 

Three  areas  differ  among  the  four  models.  They  are  the 
computational  approach,  treatment  of  focal  zones,  and  user/system 
features.  The  most  important  for  the  current  application  was  con¬ 
sidered  the  computational  approach  and  focal  zone  treatment,  since 
the  user/system  features  would  need  adaptation  from  apy  of  the 
existing  systems.  A  simplified  comparison  of  the  models  is  shown 
in  Table  2.  After  evaluation  of  the  programs,  the  TRAPS  program 
was  selected  for  two  reasons:  (a)  TRAPS  uses 


*V ir t.u’aTTy  all  BOOM  models  are  "descendants"  of  either  the  ARAP 
(Ref.  10)  or  Thomas  (Ref.  11)  models. 


a  superior  scheme  for  the  analytic  formuleition  for  ray  properties, 
(b)  the  TRAPS  program  is  the  only  program  that  allows  calculation 
of  postfocus  boom.  The  major  drawback  of  the  TRAPS  program  was 
the  lack  of  means  for  calculating  overpressures  at  focus.  It  was 
decided  that  this  could  be  remedied  by  using  the  approach  intro¬ 
duced  in  the  FOBOOM  program  for  calculations  at  focu3. 

After  decisions  were  made  to  use  the  TRAPS  program,  errors 
were  encountered  in  the  existing  TRAPS  program.  These  were 
review,  d  by  the  author  of  the  program,  Dr.  Albion  Taylor,  who  then 
developed  appropriate  corrections  to  the  program.  On  this  basis, 
the  modified  TRAPS  program  provides  answers  that  differ  from  the 
original  TRAPS  program  calculations. 
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6.0  COMPARISON  OF  SONIC  BOOM  PROGRAM  RESULTS  WITH  OTHER 
CALCULATIONS 


This  section  presents  some  comparisons  of  the  results  with  the 
modified  TRAPS  and  BOOMAP2  programs  compared  with  other  programs, 
particularly  the  FOBOOM  program,  and  with  some  field  measurements. 
Table  3  presents  some  results  for  an  F104  aircraft  in  level  flight 
and  simple  maneuvers  comparing  the  TRAPS  and  FOBOOM  programs.  The 
upper  portion  compares  results  from  the  TRAPS  and  FOBOOM  programs 
in  which  both  an  F-function  derived  from  a  nearfield  signature  was 
used  as  well  as  the  simplified  Carlson  F-function  with  Ks  =  o.07 
(see  Section  4.0)  with  field  measurements.  The  lower  portion  of 
the  table  compares  calculated  results  for  a  constant  speed  turn 
and  a  constant  acceleration  dive. 


For  the  first  case,  the  field  measurements  were  made  at 
Edwards  AFB  in  1974,  and  the  result  quoted  is  an  average  of  37 
measurements.  Using  FOBOOM,  the  effect  of  the  measured  winds  and 
non-standard  atmosphere  was  estimated  to  be  small.  The  computed 
values  shown  in  Table  3  are  for  no  winds  and  standard  atmos¬ 
phere.  The  variation  with  altitude  of  overpressure  is  shown  in 
Figure  9  and  the  signature  duration  in  Figure  10,  using  the 
F-function  from  the  near  field  signature  in  both  FOBOOM  and  TRAPS. 
The  variation  of  the  waveforms  with  altitude  is  shown  in  Figure  11 
for  TRAPS  and  Figure  12  for  FOBOOM. 


The  general  results  from  these  and  other  comparisons  indicate 
that  the  overpressure  results  calculated  by  FOBOOM  and  TRAPS  agree 
within  approximately  +5  percent.  However,  the  durations  calcu- 

1A  i.  A  U  U  .  w  jU  U  A  W  A  A  A  UM  A  a  .A  U.  A.  M  A  ^  1  .  •  A  M  A.  .HAM  A  «  MM  i.  A  M  A*  1  .  .  W  ^  f—m  1.  MM  A  U  A  fM 
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those  of  the  FOBOOM  program  by  approximately  20  percent.  The 
detailed  reasons  for  this  difference  were  not  identified  in  this 
study . 


For  the  case3  shown  in  Table  3,  there  is  essentially  no 
difference  between  the  results  from  modified  TRAPS  and  BOOMAP2. 


TAILS  ) 

CON  AAAI  SOM  OF  SOMIC  SOON  PROGRAM 
CALCULATION*  FOR  FIS')  AIRCRAFT  ID 
LSV8L  FLIC8T  AMS  SIH9LX  MAME0VER** 


Maximum 

Boom 

Flight  Condition  I 

r- 

Ground 

Overpressure ,  ref 

tXirat  ion 

Mich 

No  Altftt 

Function 

Proa  ram 

Keiaht ,ft 

+ 

- 

mate. 

I. 32 

31,200(1) 

_ 

Field  Haas 

2,330 

0.524 

... 

97 

Naar  field 

poboom 

2,300 

0.571 

-0.446 

71.2 

Naar  field 

TRAPS 

2,300 

0.556 

•  0.471 

76.8 

Carlaon 

FOBOOM 

2,300 

0.626 

-0.626 

67.9 

Carlson 

TRAPS 

2,300 

0.616 

*0.616 

93.1 

Carlson 

Simpl.  Method  2,300 

0.634 

— 

1.15 

30.000(1) 

Naar  field 

PC  BOOM 

0 

0.583 

-0.435 

91.2 

Naar  field 

TRAPS 

0 

0.SR7 

-0.426 

96.8 

Carlson 

FOBOOM 

0 

0.4  >. 

-0.660 

117.3 

Carlson 

TRAPS 

0 

0.673 

-0.673 

110.4 

1.7 

36, 000(2), (5) 

«  -  *30  * 

Carlson 

FOBOOM 

0 

0.323 

72.4 

TRAPS 

0 

0.370 

84.9 

♦  •  0* 

Carlaon 

FOBOOM 

0 

0.56) 

79.5 

TRAPS 

0 

0.559 

84.5 

a  -  .30’ 

FOBOOM 

0 

1.32 

96.2 

|  Carlaon 

TRAPS 

0 

1.24 

94.6 

♦  -  -40* 

Carlson 

FOBOOM 

0 

—  (3) 

— 

TRAPS 

0 

0.616(4) 

102.5 

1.15 

25 , 000 i  6 ) 

♦  -  0 

Carlson 

FOBOOM 

20,000 

2.07 

63.8 

15.000 

1.49 

72.0 

10,000 

1.30 

77.4 

5,000 

’.24 

81.3 

0 

1.26 

84.6 

♦  »  0 

Carlson 

TRAPS 

23,000 

2.12 

74.5 

15,000 

1.51 

82.6 

10,000 

1.30 

B7.B 

5,000 

1.20 

91.6 

0 

1.16 

94.6 

♦  -  0 

Carlson 

FOBOOM 

20,000 

2.03 

64.2 

15,000 

1.45 

72.  5 

10,000 

1.27 

78.0 

5,000 

1.22 

82.0 

0 

1.23 

05.3 

♦  -  0 

Csrlson 

TRAPS 

20,000 

2.06 

84.9 

15,000 

1.48 

83.1 

10,000 

1.27 

88.4 

5,000 

1.17 

92.  2 

n 

1.13 

95.23 

*  All  calculations  aaauia*  no  reflection  at  the 

ground 

(1) 

straight  and  lavjl  flight 

(21 

Aircraft  turning 

at  rat*  of  1. 

9425  degreee/eec 

(3) 

FOBOOM  calculates 

focus  at  701$  ft  above  ground 

(4) 

TRAPS  calculate* 

focus  at  7254 

ft  about  ground 

(5) 

«  la  tha  exumuth 

angle  of  the 

merging  rsv. 

♦  «  0  Is  dlrsctlv  below 

tha  aircraft 

(6) 

Aircraft  accalirating  at  *J  £t/aac*  at  diva 

engl*  of  27* 

58- 


FIGURE  9.  VARIATION  OF  OVERPRESSURE  WITH  ALTITUDE 
FOR  FI 04  IN  LEVEL  FLIGHT. 


FIGURE  12.  VARIATION  OF  SIGNATURE  WITH  ALTITUDE 

FOR  FI 04  IN  LEVEL  FLIGHT  -FOBOOM  PROGRAM. 
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The  effect  of  using  displacements  in  TRAPS  and  velocities  in 
B00MAP2  as  primary  flight  data  input  are  negligible  for  steady 
state  motion  of  the  aircraft. 


TRAPS,  does  not  calculate  the  overpressure  at  or  near  a  focus, 
and  for  this,  comparisons  must  be  made  between  FOBOOM  and 
B00MAP2.  Table  4  compares  FOBOOM  and  B00MAP2  results  for  an  F104 
in  level  flight,  M  *  1.154  at  25,000  ft,  accelerating  at  4 
ft/sec24  The  ground  is  at  1000  ft,  ground  reflections  are 
included  and  Carlson's  F-function  is  used.  FOBOOM  does  not 
predict  overpressure  below  a  focus,  but  the  results  in  the  focal 
region  are  within  6%,  with  the  TRAPS  results  consistently  higher 
than  those  from  FOBOOM.  The  focus  altitude  is  also  shown  in  Table 
4,  TRAPS  predicts  a  focus  approximately  80  ft  higher  than  FOBOOM. 


Table  5  compares  the  results  for  an  FI 04  in  a  level  turn  at  M 
=  1.7,  with  a  turn  rate  of  1.9425  degrees  per  second  (turn  load 
factor  *  2g),  with  the  ground  at  either  0  ft  or  1000  ft.  In  this 
case,  the  caustic  surface  is  very  steep,  with  a  focus  location  at 
1150  ft  for$*  -34.2°  and  a  slope  of  approximately  -1200  ft/degree. 


The  calculation  of  the  caustic  surface  curvature  in  the 
direction  of  the  ray  tube  collapse  is  not  always  possible  for 
every  angle  selected,  because  either  the  caustic  is  cusped  (not 
included  in  BOOMAP2)  or  the  focus  on  one  of  the  auxiliary  ray 
tubes  does  not  exist.  This  may  lead  to  gaps  in  the  results.  For 
instance,  in  Table  5  between  4  *  -32°  and  -34°,  no  overpressures 
wore  calculated  at  0  ft,  because  of  the  steepness  of  the  caustic 


surface.  However,  sufficient  rays  are  traced  to  give  a  good 
representation  of  the  focurj. 


The  BOOMAP2  and  FOBOOM  results  for  the  ground  at  1000  ft,  agree 
within  7  %  near  the  focus.  The  comparison  for  the  ground  at 
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COMPARISON  OF  FOBOOM  AND  BQOMAP2  PROGRAM  CALCULATIONS 
FOR  F104  AIRCRAFT  IN  LEVEL  FLIGHT  ACCELERATION 


Ray  Bank 
Anal© 


0° 

4° 

8° 

10° 

11° 

12° 

1 5  0 
18° 
21° 
24° 
27" 
30° 
33° 
35.3° 


Overpressure 
at  1000  ft  (psf) 

Focus  Altitude  (ft) 

BOOMAP2 

FOBOOM 

BOOMAP2 

FOBOOM 

7,448 

7.026 

985 

907 

7.040 

6.925 

1048 

964 

6.665 

6.400 

1234 

1152 

6.366 

6.186 

1375 

1296 

6.225 

5.912 

1547 

1469 

6.115 

5.912 

1547 

1469 

5.652 

5.542 

1865 

1786 

2257 

4.645 

3.622 

— 

2725 

3270 

— 

2. 845 

2.226 

— 

3894 

— 

1.706 

— 

4601 

— 

1.236 

— 

5394 

— 

0.  734 

— 

6053 

- 

NOTES : 


F104,  M  =  1.154  at  25,000  ft,  level  flight,  Acceleration 
Ground  =  1000  ft,  with  reflection  factor  of  1.0. 
Carlson's  simplified  F-function  used. 

4>  =  0  lies  vertically  below  the  aircraft. 


=  4  ft./sec2. 


Is 


r 


>*VK-V»*V  v  W-nlV*  Mr J*  IV  t  j'- '.Vtf,7i'.1VVih.V 


TABLE  5 

COMPARISON  OF  FOBOOM  AND  BOOMAP2  PROGRAM  CALCULATIONS 
FOR  F104  AIRCRAFT  IN  A  LEVEL  TURN 


Ray  Bank 


Overpressure 
at  0  ft  (ps f ) 


lie _ 

BOOMAP2 

FOBOOM 

50° 

'  0.379 

0.302 

40° 

0.594 

0.486 

30° 

0.752 

0.628 

20° 

0.890 

0.768 

10° 

1.015 

0.936 

0° 

1.137 

1.096 

-10° 

1.277 

1.248 

-20° 

1.508 

1.386 

-30° 

2.509 

2.014 

-32° 

3.829 

2.072 

-33° 

— 

2.  092 

-34° 

6.833 

8.0* 

-34.2° 

5.814 

— 

-34.4° 

— 

— 

-34.6° 

— 

— 

-34.8° 

— 

-35° 

3.773 

-36° 

2.794 

-370 

2.  245 

— 

-38° 

1.879 

— 

-40° 

1.410 

— 

-45° 

.807 

-- 

“50° 

,  460 

— 

-52.5° 

.270 

— 

Overpressure 
at  1000  ft  (psf) 
BOOMAP2  T  FOBOOM 


0.393 

0.602 

0.760 

0.896 

1.019 

1.138 

1.272 

1.484 

2.277 

3.017 

3.722 

8.300 

8.293 

8.287 

8.307 

6.739 

5.775 

3.538 

2.663 

2.154 

1.557 

.862 

.494 

.271 


0.336 

0.516 

0.652 

0.808 

0.966 

1.122 

1.296 

1.548 

2.384 

3.112 

4.054 

8.910 

8.42 


MSI 


NOTES  J 


F104,  M  =  1.154  at  36,000  ft,  level  flight.  Turn  Rate  -  1.9425°/sec 
Ground  =  0  ft  and  1000  ft,  reflection  factor  of  2.0. 

Carlson's  simplified  F-function  used. 

<J>  ■  0  lies  Directly  below  the  aircraft,  with  <J>  positive  towards  the 
center  of  the  turn. 

*  Value  from  Ref.  15. 
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0  feet  is  incomplete  in  the  focal  region.  FOBOOM  gives  a  slightly 
different  focu3  location  than  BOOMAP2,  and  hence,  the  angle  at 
which  the  focus  occurs  on  the  ground  differs  from  those  selected 
for  BOOMAP2.  As  Table  5  shows,  the  maximum  value  predicted  by 
FOBOOM  (Ref.  15)  is  8  psf,  17  percent  higher  than  BOOMAP2 
results. 
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7.0  RECOMMENDATIONS  FOR  FUTURE  WORK 


The  basic  documentation  for  the  original  TRAPS  program  (Ref. 

3)  provides  a  good  overall  discussion  of  the  major  features  of  the 
analysis  approach  incorporated  in  the  TRAPS  computer  program. 
However,  the  documentation  is  incomplete  in  that  the  analytical 
expressions  are  not  fully  described  which  makes  it  exceedingly 
difficult  to  reconstruct  the  analytic  framework  from  the  computer 
program  alone.  The  absence  of  a  detailed  mathematical  description 
handicaps  attempts  to  check  the  validity  of  the  theoretical  model 
and  to  extend  the  analytic  approach  to  cover  features  missing  from 
the  original  program.  This  lack  is  particularly  unfortunate, 
because  the  basic  TRAPS  analytical  approach  gives  every  evidence 
of  being  superior  to  those  utilized  in  earlier  ray  tracing 
models. 

It  is  therefore  recommended  that  the  theoretical  basis  for  the 
TRAPS  program  be  developed  and  documented  in  detail.  The  best 
person  to  do  this  would  be,  of  course,  the  original  author  of  the 
TRAPS  program,  Dr.  Albion  Taylor.  Such  documentation  would  not 
only  provide  a  basis  for  the  detailed  evaluation  of  the  TRAPS  pro¬ 
gram  in  handling  various  types  of  sonic  boom  situations,  but  would 
also  provide  a  basis  for  possible  extension  of  the  TRAPS  program 
to  cover  sonic  boom  situations  that  were  not  of  immediate  interest 
in  the  current  program,  but  which  may  be  of  vital  importance  in 
other  applications.  For  example,  the  capabilities  of  the  TRAPS 
program  to  handle  rays  that  have  risen  to  a  high  altitude  and 
returned  to  the  ground  were  not  of  interest  in  the  current  appli= 
cation,  but  may  be  vital  in  future  applications  of  the  program  of 
interest  to  the  Air  Force  as  well  as  to  other  users.  Also,  the 
lack  of  analytic  formulation  limits  the  rigorous  comparison  of  the 
program  predictions  with  the  results  of  other  programs  or  of  field 
experiments . 


7 . 2  Extension  of  the  TRAPS  Analytic  Framework  to  Include  Focus 
Slqnaturaa  and  Overpressures 

None  of  the  existing  ray  theory  programs  (including  the  TRAPS 
program)  predicts  the  overpressure  and  wave  signature  at  focus  as 
part  of  their  original  mathematical  development.  The  available 
programs  that  can  handle  this  situation  (FOBOOM  and  the  modified 
TRAPS  program)  incorporate  an  approach  to  estimating  the  signature 
and  overpressure  that  is,  more  or  less,  "grafted"  on  the  basic 
program.  This  grafting  process  and  the  assumptions  inherent  in  it 
lead  to  potentials  for  inexact  calculations  and,  perhaps, 
oversimplifications.  It  is  recommended  that  serious  effort  be 
given  to  extending  the  theoretical  basis  of  the  original  TRAPS 
program  to  include  a  consistent  analytic  model  for  computing  the 
overpressure  and  wave  signatures  at  sonic  boom  focu3  locations. 
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APPENDIX  A 

Sample  BOO MAP 2  Program  Output 

This  appendix  provides  a  sample  of  the  300MAP2  computer  program 
output  based  upon  analysis  of  the  MOAOPS  library  tapes  for  Luke 
Air  Force  Base.  Figure  A-l  lists  the  missions  included  in  the 
analysis.  At  the  end  of  the  figure,  the  number  of  supersonic 
flights  is  listed,  together  with  the  number  of  boom-producing 
flights.  These  figures  are  based  upon  the  telemetered  Mach  num¬ 
bers,  applying  criteria  derived  from  the  steady  flight  Carlson 
equations. 

At  the  end  of  Figure  A-l,  the  total  supersonic  time  and  total 
boom-producing  time  are  also  listet  One  set  of  figures  is 
based  upon  the  telemetered  Mach  number.  The  other  set  of  fig¬ 
ures  is  based  upon  a  Mach  number  calculated  from  the  ground 
velocity  and  the  standard  day  temperature.  Both  set  of  figures 
are  based  upon  criteria  assuming  steady  flight  ( Carlson  equa¬ 
tions).  All  of  the  statistics  for  boom-producing  flights  and 
times  listed  at  the  end  of  Figure  A-l  are  approximate  since  they 
are  based  upon  steady  flight  assumptions.  In  the  detailed 
BQ0MAP2  analysis,  using  the  modified  TRAPS  ray-tracing  program, 
the  total  boom-producing  time  may  be  quite  different  since  the 
B00MAP2  calculations  do  take  into  accoui  t  accelerations  and 
turns . 


Figure  A-2  presents  various  statistics  for  the  geographic  loca¬ 
tion,  aircraft  altitudes,  Mach  numbers  and  sound  levels  (with 
sound  levels  calculated  from  the  Carlson  equations).  Figure  A-3 
presents  information  on  the  geographic  location  of  supersonic 
flights . 


A-l 


Figure  A-4  provides  a  sample  of  the  "scratchpad"  for  one  flight 
that  generated  a  focus  sonic  boom.  Because  of  the  crude  con¬ 
touring  program  employed/  the  resulting  display  of  the  overpres¬ 
sure  contours  is  simplified. 

Note  that  the  geographic  location  of  the  calculated  maximum 
overpressure  is  shown  on  the  scratchpad  plot.  The  information 
on  the  geographic  location  of  the  maximum  overpressure  is 
collected  froia  each  scratchpad  and  used  to  generate  a  map 
showing  the  location  of  all  focus  booms  for  the  MOAOPS  library 
set  under  study  (see  Figure  A-7). 

Figure  A- 5  shows  the  geographic  location  of  the  flight  tracks 
for  all  supersonic  aircraft  activity  (as  based  upon  reported 
Mach  number).  Figure  A-6  shows  the  flight  tracks  that  were 
likely  to  produce  sonic  booms  that  reached  the  ground  (based  on 
the  cut-off  equations  of  Carlson). 

Figure  A-7  shows  the  geographic  location  of  the  maximum  over¬ 
pressures  resulting  from  focused  sonic  booms.  For  each  of  the 
geographic  locations  shown,  there  is  a  corresponding  scratchpad 
which  provides  information  about  the  flight  that  produced  the 
focus  boom. 

Figure  A-8  shows  sonic  boom  noise  contours  for 'Luke  AFB .  Figure 
A-8  shows  the  average  C-weighted  sound  exposure  level  contours. 
Contours  of  average  maximum  overpressure  in  pounds  per  square 
foot  (psf)  can  also  be  produced  by  B00MAP2. 


SOURCE  LISTING 


TITLE  LUXE  HUE.  Cf  ENTIRE  LIBRARY 
....  SITE  LUXE 

3:  DATE  ALL 

A:  TINE  ALL 

5:  AIRCRAFT  ALL 
S;  KACHTRK  W50DDO 
/!  BOCNTRK  W'f.OOO 

H;  CONTOUR  r  .IL, 440000  ,  93.,  TOO.,  103. p  110.,  IIS.,  120. 

T:  STATS 


1/  1 


i..a.P .  c  :...»,alvaili»»„iifiKMiiii>»i,raaaaainaK,iiUi,i,u,iT  .........a 

I  SITE  I  l.\  '  “SE  I  OATE  I  TIM!  IAIRCSAFTI  TAIL  I  •  I 

I  LOCATION  I  I  nrYNHDO-TYMKOD]  !  CHHNH-HHW]  I  TYPE  I  NUMBER  !  I 


I . . . . . . I 

I  LUKE  I  ALL  I  ALL  I  ALL  I  ALL  t  ALL  !  1  I 


..a-jraair  '.-.^'.nraaaaaascanxxaaasaxx.iaKaaaaKaiaxitxaaaaaaaaaxKaiMiaaaaaaaamaiawnBaa 
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4  FIGURE  A~1 .  CONTINUED 


NUMBER  Of  BOOM  PRODUCING  SORT  I  ESC  FLIGHTS):  10 

USING  MACH  NO  CALCULATED  FROM  GROUND VELOCITIES 
TOTAL  SUPERSONIC  TIME  ■  1528  SECONDS 

TOTAL  BOOM  PRODUCING  TIME  -  393  SECONDS 

USING  TELEMETERED  MACH  HO  CALCULATED  FROM  AIRSPEED 
TOTAL  SUPERSONIC  TIME  -  2955  SECONDS 

TOTAL  BOOM  PRODUCING  TIME  «  874  SECCNOS 
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rLIGTfT  TRACK  SEGMENTS  OF  SUPERSONIC  AIRCRAFT 
ACTIVITY  (MACH  »  1) 
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FIGURE  A- 5.  SUPERSONIC  AIRCRAFT  FLIGHT  TRACK  MAP  FOR 
LUKE  AFB 
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I  DXATION  OF  FOCUS  BOCM  MAXIMUM  OVERPRESSURES 

SCALE:  1  INCH  »  40000  FEET 
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FIGURE  A- 7.  LOCATION  OF  FOCUS  SONIC  ROOM  MAXIMUM  OVERPRESSURES 
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FIGURE  A-8.  AVERAGE  C -WEIGHTED  SOUND  EXPOSURE  LEVELS  FOR 
I.UKE  AFB  LIBRARY 


